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Quantum chaos and chiral symmetry at the QCD and QED phase 
transition 
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We investigate the eigenvalue spectrum of the staggered Dirac matrix in SU(3) gauge theory and in full QCD 
as well as in quenched U(l) theory. As a measure of the fluctuation properties of the eigenvalues, we consider 
the nearest-neighbor spacing distribution. We find that in all regions of their phase diagrams, compact lattice 
gauge theories have bulk spectral correlations given by random matrix theory, which is an indication for quantum 
chaos. In the confinement phase, the low-lying Dirac spectrum of these quantum field theories is well described 
by random matrix theory, exhibiting universal behavior. Related results for gauge theories with minimal coupling 
are now discussed also in the chirally symmetric phase. 
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1. Quantum chaos 

The eigenvalues of the Dirac operator are of 
great interest for the universality of important 
features of QCD and QED. On the one hand, 
the accumulation of small eigenvalues is, via the 
Banks-Casher formula ||l[, related to the sponta- 
neous breaking of chiral symmetry. On the other 
hand, the fluctuation properties of the eigenvalues 
in the bulk of the spectrum can be described by 
random matrix theory (RMT), see Ref. 0. For 
example, the nearest-neighbor spacing distribu- 
tion P(s), i.e., the distribution of spacings s be- 
tween adjacent eigenvalues on the unfolded scale, 
agrees with the Wigncr surmise of RMT. Ac- 
cording to the Bohigas-Giannoni-Schmit conjec- 
ture P], quantum systems whose classical coun- 
terparts are chaotic 0] have a nearest-neighbor 
spacing distribution given by RMT whereas sys- 
tems whose classical counterparts are integrable 
obey a Poisson distribution, Pp(s) — e^^. There- 
fore, the specific form of P{s) is often taken as a 
criterion for the presence or absence of "quantum 
chaos" . 

In RMT, one has to distinguish several univer- 
sality classes which are determined by the symme- 
tries of the system. For the case of the QCD Dirac 
operator, this classification was done in Ref. [g|. 
Depending on the number of colors and the rep- 
resentation of the quarks, the Dirac operator is 
described by one of the three chiral ensembles of 



RMT. As far as the fluctuation properties in the 
bulk of the spectrum are concerned, the predic- 
tions of the chiral ensembles are identical to those 
of the ordinary ensembles [g|. In Ref. [||, the 
Dirac matrix was studied for color-SU(2) using 
both staggered and Wilson fermions which corre- 
spond to the chiral symplectic (chSE) and orthog- 
onal (chOE) ensemble, respectively. Here 0, we 
study SU(3) with staggered fermions which cor- 
responds to the chiral unitary ensemble (chUE). 
The RMT result for the nearest-neighbor spacing 
distribution can be expressed in terms of so-called 
prolate spheroidal functions, see Ref. ||]. A very 
good approximation to P{s) is provided by the 
Wigner surmise for the unitary ensemble. 
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We generated gauge fleld conflgurations using 
the standard Wilson plaquette action for SU(3) 
with and without dynamical fermions in the 
Kogut-Susskind prescription. We have worked 
on a 6^ X 4 lattice with various values of the in- 
verse gauge coupling /3 — G/g^ both in the con- 
finement and deconfinement phase. We typically 
produced 10 independent equilibrium configura- 
tions for each /3. 

The Dirac operator, p = ^ + ig^, is anti- 
Hermitian so that the eigenvalues A„ of ip are 
real. Because of {p, 75} = the nonzero A„ occur 
in pairs of opposite sign. All spectra were checked 





Figure 1. Nearest-neighbor spacing distribution P{s) for full QCD on a 6'^ x 4 lattice in the confinement 
phase (left) and in the deconfinement phase (right) compared with the random matrix result (solid lines). 
There are no changes in P{s) across the deconfinement phase transition. 
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Figure 2. Nearest-neighbor spacing distribution P{s) for U(l) gauge theory on an 8'^ x 6 lattice in the 
confined phase (left) and in the Coulomb phase (right). The theoretical curves are the chUE result, 
Eq. (j^), and the Poisson distribution, Pp(s) = e"**. 



against the analytical sum rules ^^^ A„ = and 
Sa >o ^n — 3V^, where V is the lattice volume. 
To construct the nearest-neighbor spacing distri- 
bution from the eigenvalues, one first has to "un- 
fold" the spectra g. 

Figure compares P{s) of full QCD with Nf = 
3 flavors and quark mass ma = 0.05 to the RMT 
result. In the confinement {(3 — 5.2) as well as 
in the deconfinement (/3 = 5.4) phase we ob- 
serve agreement with RMT up to very high /3 
(not shown). The observation that P{s) is not 
infiuenced by the presence of dynamical quarks is 
expected from the results of Rcf . |g] , which apply 
to the case of massless quarks. Our results, and 
those of Ref . , indicate that massive dynamical 
quarks do not affect P{s) either. 

We have also investigated the staggered Dirac 



spectrum of 4d U(l) gauge theory which corre- 
sponds to the chUE of RMT. At /3c « 1-01 U(l) 
gauge theory undergoes a phase transition be- 
tween a confinement phase with mass gap and 
monopole excitations for (3 < (3c and the Coulomb 
phase which exhibits a massless photon for [3 > 
(3c §■ As for SU(2) and SU(3) gauge groups, 
we expect the confined phase to be described 
by RMT, whereas free fermions yield the Pois- 
son distribution. The question arose whether the 
Coulomb phase would be described by RMT or 
by the Poisson distribution |1^. The nearest- 
neighbor spacing distributions for an 8'^ x 6 lat- 
tice at /3 — 0.9 (confined phase) and at (3 — I.l 
(Coulomb phase), averaged over 20 independent 
configurations, are depicted in Fig. |[ Both are 
wefl described by the chUE of RMT. 



No signs for a transition to Poisson regular- 
ity are found. The deconfinement phase tran- 
sition does not seem to coincide with a transi- 
tion in the spacing distribution. For very large 
values of (3 far into the deconfinement region, 
the eigenvalues start to approach the degener- 
ate eigenvalues of the free theory, given by A^ = 
^ ;^sin^(27rn^/L^)/a^, where a is the lattice 
constant, L^ is the number of lattice sites in the 
//-direction, and n^ = 0, ...,L^ — 1. In this 
case, the spacing distribution is neither Wigner 
nor Poisson. It is possible to lift the degeneracies 
of the free eigenvalues using an asymmetric lat- 
tice where L^;, Ly, etc. are relative primes and, 
for large lattices, the distribution is then Poisson, 
Pp(s) = e-", see Fig. |. 
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Figure 3. Nearest-neighbor spacing distribution 
P{s) for the free Dirac operator on a 53 x 47 x 43 x 
41 lattice compared with a Poisson distribution. 

2. Chiral symmetry 

We have continued our investigations with a 
study of the distribution of the small eigenval- 
ues in the whole phase diagram. The Banks- 
Casher formula Q relates the Dirac eigenvalue 
density p{X) at A = to the chiral conden- 
sate, E = \{ip^)\ — limg^o limy^oo 7rp(e)/V'. 
The microscopic spectral density, Ps(2) = 
limy_^oo P {.zjy^ /^S, should be given by the 
appropriate prediction of RMT ||ll|, which also 
generates the Leutwyler-Smilga sum rules Q. 

We present results in Fig. ^ for SU(3) theory 
and the staggered Dirac operator on a 4^ lattice 
from 5000 configurations for /3 = 5.4 and 3000 
configurations for /3 = 5.6 and /3 = 5.8 around 




the critical temperature /3c ~ 5.7. In the confine- 
ment phase, both the microscopic spectral density 
Pi(z) and the distribution P(Amin) of the small- 
est eigenvalue agree with the RMT predictions of 
the chUE for topological charge z^ = Q . 

Our analog presentation for U(l) theory is from 
10000 configurations on a 4* lattice around the 
critical coupling. In the left plot of Fig. g a com- 
parison with RMT for the microscopic spectral 
density Ps(z) yields again quite satisfactory agree- 
ment in the confinement. The analytical RMT re- 
sult for the (quenched) chUE and i/ = is given 
by p,{z) = z [J2(z) + Jf (z)]/2, where J de- 
notes the Bessel function. The chiral condensate 
E can be obtained by extrapolating the histogram 
for p(A) to A = and using the Banks-Casher 
relation |14|. The right plot in Fig. exhibits 
the distribution P(Amin) of the smallest eigen- 
value, being in the chirally broken phase in ac- 
cordance with the prediction of the (quenched) 
chUE of RMT for topological charge v — Q, 

P{Xn.in) = (l/E)2(A„,in/2) exp(- (FEAn^in) V4) • 

The quasi-zero modes which are responsible for 
the chiral condensate E ^ build up when we 
cross from the deconfinement into the confined 
phase. Figures |j and [^ demonstrate that both 
p(A) and P{Xmin) plotted with varying /? on iden- 
tical scales, respectively, can serve as an indicator 
for the phase transition. 

In Fig. n we turn to a discussion of the spec- 
trum in the quark-gluon plasma and Coulomb 
phase. From RMT a functional form of p(A) = 
C(A - Ao)^™"*"^/^ is expected at the onset of the 
eigenvalue density [|5[. A fit to the data in the 
regime up to A = 1 yields m — 0.081(1) for the 
non-Abelian and ni — 0.084(1) for the Abelian 
theory, in agreement with recent studies |16[ . 
This suggests that both theories correspond to 
universality class to = 0. For this class a micro- 
scopic level density involving the Airy function 
can be deduced from RMT ||l^. A rescaling of 
our data from the 4* lattice to this functional 
form is not satisfactory for both theories |16| . 
The reason might be (i) that the lowest eigen- 
value is still influenced by quasi-zero modes from 
the confinement and (ii) that already the first 
eigenvalues lie above an analogue of the Thou- 
less energy. With increasing lattice size these ef- 
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Figure 4. Density p{X) of small eigenvalues (left) and distribution P(Xmin) (right) for SU(3) gauge theory 
on a 4^ lattice across the phase transition. 
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Figure 5. Density p{X) of small eigenvalues (left) and distribution P{Xmin) (right) for U(l) gauge theory 
on a 4 lattice across the phase transition. 
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Figure 6. Fit of the spectral density to p(A) = C(A - Ao)2"+^/2 in SU(3) at /3 = 5.8 (left) and in U(l) at 
/3 — 1.10 (right). The contribution of the smallest eigenvalue, the 11^ eigenvalue and the 21®^ eigenvalue 
is inserted. 



fects should decrease. Further, we checked that 
the high-end of the spectrum behaves similar to 
the low-end and thus deviates from the results of 
the ordinary UE of RMT [0 . Again the ques- 
tion arises whether larger lattices possess spectral 
edges closer to the microscopic semi-circle univer- 
sality. 

3. Conclusions 

The aim of this contribution was to work out 
two different types of universalities inherent in 
quantum field theories with a covariant deriva- 
tive and their interpretation in terms of RMT. 
The first type concerns the bulk of the spec- 
trum of the Dirac operator. The nearest-neighbor 
spacing distribution P{s) agrees with the RMT 
result in both the confinement and the decon- 
finement phase of pure gauge theory and of full 
QCD, except for extremely large values of (3 
where the eigenvalues are known analytically. 
The nearest-neighbor spacing distribution of 4d 
U(l) quenched lattice gauge theory is described 
by the chUE of RMT in both the confinement and 
the Coulomb phase. Even in the deconfinement 
phase, gauge fields retain a considerable degree of 
randomness, which apparently gives rise to quan- 
tum chaos in these theories. 

The second type of universality concerns the 
low-lying spectra of the Dirac operators of both 
QCD and QED. In all cases considered, one finds 
that in the phase in which chiral symmetry is 
spontaneously broken the distribution P(Ainin) 
and the microscopic spectral density Ps{z) are de- 
scribed by chiral RMT. When chiral symmetry is 
restored one has to rely on ordinary RMT with 
the related space-time symmetries, but one finds 
for our lattice size only universal behavior of the 
macroscopic density /9(A). 
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